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ABSTRACT 



This thesis is comprised of four independent computer 
programs and their related operating instructions. Each of 
these programs focuses on a particular facet of tactical 
missile design. The topics covered include guidance and 
trajectory calculations, rocket motor propulsion sizing, 
warhead design, and aerodynamic coefficient determination. 
The programs are developed from accepted mathematical proce- 
dures and are tailored to optimize operator interaction for 
educational purposes. This thesis is intended to be uti- 
lized as a supplement to the thesis Tactical Missile Concep - 
tual Design by D.R.Redmon, Naval Postgraduate School, Sep- 
tember 1980. 
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I. INTRO DOCTIQN 



The programs contained in this thesis were assembled 
expressly to supplement the work of Dan Redmon in his the- 
sis, Tactical Missile Conceptual Design . Two of the 
programs, LPATH, the trajectory model, and LAER01 , the aero- 
dynamic coefficients program, originated in Redmon' s thesis 
and were expanded/modified for use on the Naval Postgraduate 
School's new IBM 370 computer system. The other two pro- 
grams were generated for this thesis and utilize the proce- 
dures and principles outlined by Redmon. 

The specific intention of these programs is to provide 
students of tactical missile engineering and design with a 
method of solving complex mathematical routines rapidly and 
interactively. Each of the programs request data which are 
likely to be used as design parameters for the topics con- 
cerned. The programs also allow repeated operation with 
input alteration capability, allowing the user to personally 
optimize his design. This approach was chosen to allow stu- 
dents to understand the relationships various input parame- 
ters have with the final solutions. 
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II. TRAJECTORY MODELS 



A. DESCRIPTION AND ORIGIN OF PROGRAM 

This program applies the principles of missile guidance 
laws to the terminal phase (the last 5 to 10 seconds) of a 
missile trajectory in order to determine the maximum normal 
acceleration on the missile for a given scenario. Of the 
three general guidance law categories, pursuit guidance is 
not included in the program capability. It has been found 
that pursuit guidance invariably produces a tail-chase situ- 
ation, greatly reducing an anti-air missile's effectiveness 
against maneuvering targets of similar speed characteris- 
tics. Line-of-sight guidance and proportional navigation 
guidance are both options of the program. 

Figure (II— 1 ) shows a typical encounter geometry as 
required for this program. The encounter is considered to 
occur entirely within a two- dimensional plane. No differ- 
entiation is required or assumed concerning the orientation 
of the encounter plane. The plane may be at any angle to 
the horizontal as desired by the program user. The refer- 
ence direction is an arbitrary choice by the program user. 
The angles shown are positive in value, however the program 
does not reguire positive angles. If TAL were 150 degrees, 
it could also be entered as -210 degrees. The IRA term 
represents the initial range to the target. 

Tangential velocities (air speeds) of rhe missile and 
the target are considered by the program to be constant 
throughout the problem. Since the program concerns itself 
with only the final moments of a trajectory, this is a rea- 
sonable consideration. Target normal accelerations, when 
specified by the user, are also held constant throughout rhe 
problem for the same reason. The missile normal 
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Program Variables 
j3=L0S (line of sight angle) 
target alpha) 
missile alpha) 
target speed) 
missile speed) 
initial range to 
target from missile 
launch point) 
(initial range to 
missile from 
launch point) 



<*t=TAL 
atm=M AL 
7t=TSP 
Vm=M SP 
Rt=IRT 



Rm=IR« 



/ 

/ 




/ 




Reference direction 



Figure (II-1) . Encounter geometry 
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acceleration, however, is a function of the guidance laws 
and will vary according to the resulting flight path 
requirements. A constraint on the encounter is that the 
missile normal acceleration must be zero at the start of the 
problem. 

The program "flies*' both the missile and the target as 
simple points in space with no consideration given to aero- 
dynamics. The missile will always strike the target dead 
center when it is given the proper speed advantage for the 
encounter since there is no provision for statistical miss 
analysis. The unbounded nature of the missile normal accel- 
eration allows the missile to turn as rapidly as necessary 
to hit the target. 

This program analyzes the given encounter by time incre- 
ment calculations. As is the case with all integrations 
conducted by incremental steps, the accuracy of the results 
is a function of the increment size. The results will 
become increasingly accurate as the time increment is made 
smaller. However, as the time increment decreases in size, 
the length of the output becomes increasingly longer. The 
user must balance the desire for accuracy against the amount 
of time he wishes to spend on the computer terminal. 

The primary output is a tabular listing of the missile 
and target coordinates at each time increment. The coordi- 
nate frame is established within the encounter plane with 
the abscissa oriented along the reference direction. The 
problem stops once the missile has passed its closest point 
of approach to the target. Output then provides the time to 
intercept from time of problem initiation and the maximum 
acceleration the missile was required to endure. A Versatec 
plot of overlaying successive encounters is an optional 
output. 
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This program originated as two seperate BASIC programs 
written by Eedmon (Ref. 1] for use on the HP 9830 desktop 
calculator. It was subsequently translated into FORTRAN IV 
for use on the IBM 370 computer system. 



B. USER INSTRUCTIONS 



If it is desired to abort the operation of this program 
prematurely, two methods are available. When the program is 
waiting for data entry, press (ENTER[ . When the program is 
not waiting for data entry but is processing, type "HX" and 
press j ENTER | . Both actions will return the terminal to CMS 
mode. 

When the screen becomes full, or "MORE...." appears in 
the status area, clear the screen by pressing ALT and 
) CLEAR ( simultaneously. At several points in the program, 
the user will be directed to "CLEAR SCREEN AND INPUT 0". 

This is to provide proper positioning of the output on the 
screen for ease of reading. If any other symbol than 
“EF®” should appear in the lower left of the screen, 
press j RESETj and continue. r-i 

1. Turn the terminal on with the red lol switch. 

2. When the large "NPS" appears after about 30 seconds, 
press | RESET! , then press [ENT ER} . 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nnnnP", where nnnn is 
your 4-digit user number, then press 1 ENTER 1 . 

4. you will now be asked for your password. Type it in 
(the characters will not appear on the screen) , then press 

\ ENTER 1 . 



5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 
does not, type "GLOBAL TXTLI3 FORT MOD2 MOD2EEH" , then press 

1 enter! • 
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6. Type "LliJK TO xxxxP 191 AS 192 SB", where xxxx is 
the 4-digit user number for the project file, then press 
1 enteb | . 



7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 
then press {ENTEB | . 

] 



ENTEB 



8. Type M ACCESS 192 B" and press 

9. Press } ALT j and | CLEAB ) simultaneously to clear 



screen. 

10. Type "LPATH" and press | EM TEH) . 

11. Input the following data as it is reguested on the 
screen by typing the data and then pressing | ENTEB 



Ensure 



that the data is input as either decimal or integer as spe- 
cified. The terms in parenthesis below are the program 
variable names. 



Parameter 


Units 


Value ranqe 


Trajectory option (TITLE) 


none 


0=Line-of- sight 
1 =ProDortional 
navigation 
2=Both 


Time increment (DEL) 


seconds 


Larger than .0001 
times the problem 
time 


Navigation constant (NAV) 


none 


2.5 to 6.5 


LOS Angle (LOS) 


degrees 


0.0 to 360.0 


Initial target range 
from the launch site (IBT) 


meters 


Larger than the 
missile range 


Target speed (TSP) 


m/sec 


Larger than 0.0 


Target normal 
acceleration (TAC) 


m/sec/sec 


Positive is to 
target's left 


Initial target angle to 
line of sight (TAI) 


degrees 


0.0 to 360.0 the 


Missile speed (MSP) 


m/s ec 


Larger than 0.0 


Initial missile range 
from launch site (IBM) 


meters 


Such that impact 
occurs in less 
than 10 seconds 



12. If you desire to rerun the problem, or want to run a 
new problem, answer the questions appropriately when asked 
by the terminal after the execution of the current problem. 
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13. To receive the printout and plot of your encounters, 
answer "no" to rerunning or restarting the problem when 
asked by the terminal and follow the directions presented on 
the screen. 

14. Upon completion of the program, type "LOGOFF" and 



press | ENTER| . 

15. Turn the terminal off with the red 



switch. 



C. EXAMPLE PBOBLEMS 

1 . Example II-A . Line-of -sjg ht Non-maneuvering Cross - 
ing Target 




Figure (II-2) . Non-maneuvering crossing target 

Rt=10000 meters 
Bm=9000 meters 
do=30. 0 degrees 
0<to=-90.0 degrees 
at=0 

Vt=200 meters/second 
Vm=800 meters/second 
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Figure (II- 0. Versatec plot of example II- 
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Table (II- 1) is the corresponding computer output for 
this encounter. As indicated, the missile maximum normal 
acceleration is 

am= -32.00 m/sec/sec or -3.26 g's. 

Figure (II-3) is the Versatec plot of the engagement. 

2 . Example II-B . Proportional Navigation Maneuvering 
crossing Target 




Launch 

point 




Figure (II-4). Maneuvering crossing target 

Rt=10000 meters 

Hm=9000 meters 

#o=30.0 degrees 

0<to=-9 0.0 degrees 

at = 156. 8 m/sec/sec (16 g's) 

Vt=200 meters/second 
Vm=800 meters/ second 
Navigation constant=3 . 06 
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Figure (II-5) . Versatec plot of example II- 
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Table (II-2) is the computer output for this encounter. 
The missile maximum normal acceleration is 

am= 283.94 m/sec/sec or 28.95 g's. 

Figure (II-5) is the Versatec plot of the engagement. 
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D. PROCEDURAL FLOWCHART 



< 



Input data 

: i::::: 



Do you want just proportiona 
navigation? 






!> 



YES 



Mo 



Establish initial angular 
relationships and determine 
initial 

X-Y coordin tes for missile 
and target 



I 



Increment iteration counter 

I 



< Have you exceeded 9999 itera-^^N* 

tions? 

Print error message 
[ END ] 



Mo 



Calculate angular relation- 
ships between missile and tar- 
get 



Calculate normal acceleration 
on missile 



Determine range between mis- 
sile nd target 



Calculate new X-Y coordinates 
for missile and target based 
on speeds 

and angular relationships 



I 



Increment time 

YES ^ 1-_— V 

— — Is range greater than zero? 
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kIc 



<: 



t 

Print output 



Do you want just line-of- 
sight ? 



I 

[ END J 




Establish initial anaular 
relationships and determine 
initial X-Y 

coordinates far missile 5 tar- 
get 



I 



H Increment iteration counter 



< H 



I 



Have you exceeded 9999 itera- 
tions 



I 



VES 



Print error message 
END 



Calculate new X- 
for missile & tar 
soeed 

and angular' rel 



pc 

get based on 
ationships 



I 



Determine range 
sile & ta 



1 



Increment time 



I 



Calculate miss 
acceleration regu 
on chan~ 
angular rela 



I 



Calculate new an 
tix> nshi 



T 



gular rela- 
ps 



between mis- 
rget 



ile normal 
irment based 
e in 

ion ships 
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JSS <C III range greater than zero? :> 

Print output 
[ END ] 



E. PROGRAM CHANGES 

1 . Language Translating . 

The two BASIC pro rams contained in Reference 1 were 
translated into standard FORTRAN IV. 

2 . Program Condensation . 

The two individua programs were combined to form a 
single integrated routine which allows the user to choose 
either of the two guidance laws or both for a given encoun- 
ter. The two original programs were meshed such that only 
the actual guidance law algorithms are separate routines, 
all input and output rout nes are now common. 

3 . Input and Output Facility . 

The data input instructions were modified to maxim- 
ize user facility on the ISM 370 System 3278 terminals. 
Completion of data input is now followed by data feedback 
for user verification prior to actual program execution. 

Data output has been expanded to provide data deliv- 
ery to three destinations; to the user's terminal for obser- 
vation, to a printer file for hardcopy duplication of the 
terminal display, and to a plot data file for subsequent use 
by the plot program. The destinations are options chosen by 
the user for each execution of the program. Dp to nine 
different problems can be printed and plotted each time the 
program is entered. 
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The program can now be rerun multiple times without 
exiting and re-entering each time. The user has a choice of 
either rerunning the same problem or initiating a new prob- 
lem completely. 

4 . Plot Program . 

PATHPLOT FORTRAN was developed to produce a Versatec 
plot picture of the encounter. It will produce a single 
plot each time the program is entered and will plot up to 9 
engagements in an overlaying manner. This format was 
adopted to allcw comparisons of successive input data modi- 
f ications. 

5. Data Overcapacity Check . 

If the user initiates a problem requiring more than 
9999 iterations, the program will stop. The user will be 
notified of the error and given the opportunity to rerun the 
problem. 

6 . Initial Missile Acceleration . 

Initial missile acceleration was removed as a user 
input variable and established as zero. Due to the mathe- 
matical nature (i. e. , no physical constraints) of the pro- 
gram, any "wrong" accelerations of the missile in the ini- 
tial state were immediately "corrected" by the algorithm. 

The model is better served by providing no initial accelera- 
tions. 

7 . Theta Angles . 

Both the target and the missile theta angles (the 
angles between the velocity vectors and the reference direc- 
tion) were removed as user input variables. The program now 
calculates the theta angles from other input data, reducing 
redundancy and possible contradiction of input data. 

3 . Initial Conditions Perspective . 

Originally, the missile guidance calculations 
started at t=0. Specifically, the anchor point for the line 
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of sight solution was the point where the missile commenced 
the problem/ whether or not that was on the launch site. 

This produced a situation removed from reality where the LOS 
anchor point should be at the fire control location, usually 
at or near the launch site. & similar, though less pro- 
nounced, condition existed for the proportional navigation 
solution. The program was modified to provide proper posi- 
tioning of the external guidance reference. As a result of 
the modification, an additional output is the correct lead 
angle for the missile at the start of the problem. This 
angle is based on the assumption that neither the target nor 
the missile have maneuvered prior to time t=0. 
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III. WARHEAD DESIGN 



A. DESCRIPTION AND ORIGIN 

This program develops a warhead using the same methods 
as presented by D Redmon (Ref. 1], However, its capability 
is somewhat greater and it applies the relationships in a 
slightly different manner. This program was written for 
this thesis in FORTRAN IV for use on the IBM 370 computer 
system 3278 terminals. It is configured expressly for oper- 
ator-computer interaction. 

It starts with an initial input of data from which a 
table of' fragment initial velocities is generated and pre- 
sented to the user. From this, the user selects values for 
fragment mass and impact velocity and another table is gen- 
erated and displayed. From this one, the desired probabil- 
ity of hit given a detonation is selected and the final 
solution is produced. At various points during the opera- 
tion of the program, the user will have opportunities to 
alter or revise selected parameters. 

The program is limited to a cylindrical warhead with 
either a solid or hollow core. The fragments- are required 
to be square in shape and are sized by the program. Figure 
(III— 1) illustrates the shape of the warhead and the loca- 
tion of various input and output quantities. 

Initially, the target altitude is used to determine the 
atmospheric density, temperature and speed of sound. These 
values are, in turn, used to calculate required initial vel- 
ocities for the fragments. The program is preloaded with 
various values for the fragment mass and impact velocity, 
which are used to generate the table of initial velocities. 
These velocities come from the following relationship: 

V (hit) = Vi (exp (-ks) J 
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k=2^ fkkCd 

where V (hit) is the impact velocity, Vi is the initial vel- 
ocity, s is the kill radius, m is the fragment mass, is 
the atmospheric density, A is the plan area of the fragment, 
and Cd is the drag coefficient of the fragment. 

The ballistic limit velocity is calculated for the vari- 
ous presized fragments and is provided as a reference when 
choosing an appropriate impact velocity. The ballistic 
limit velocity is that velocity at which one half of the 
fragments will penetrate the target's skin and the other 
half will not. The empirical relationship, developed by A. 
E. Fuhs £Hef. 6], presents the ballistic limit velocity as a 
function of the fragment size to skin thickness ratio. His 
function dealt with steel fragments impacting an aluminum 
plate. His results were modified to qualitatively account 
for different skin materials and fragment densities. 

Next computed is the fragment spray angle and the criti- 
cal miss distance. The spray angle is a function of the 
initial velocity, the detonation velocity and the warhead 
length. The critical miss distance is defined as the range 
where the fragment spray exactly covers the entire target. 
The critical miss distance is used by the program to sepa- 
rate the two functions which determine the average number of 
hits received by the target. The program assumes the target 
is always centered within the fragment spray and aligned 
perpendicular to the central ray of the spray. 

A selection of warheads is then produced, one for each 
of a preloaded set of Pd's (probability of a hit given a 
detonation) to provide the user with a Pd versus warhead 
weight/size trade off comparison. This sizing process is 
based entirely on the following relationship: 

Pd=1 -exp (-a) 
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where a is the average number of hits. The variable a, as 
shewn by Redmon (Ref. 1] , is a function of the cube of the 
warhead radius. The user then chooses a desired Pd which, 
in turn, produces the final warhead sizing. 

Values that are shown in tabular form for user selection 
and input into the program are not limited to those dis- 
played. Any value in between the displayed values or com- 
pletely outside of the value range may be used. The one 
exception is that Pd can never be selected to be larger than 
.999 and may even, if forced by the program limitations, be 
reguired to be lower. Since ultimately in this algorithm, 

Pd is a function of the warhead radius, the maximum value 
for Pd may be reduced in order to keep the radius within the 
original input parameters. The user is notified if this 
condition occurs. 

Useful reference information for some common explosives 
and case materials is contained in the following tables. 



Table (III-1). Characteristics of common explosives 

EX |i sl ” 



RDX 


. 0596 


9300 . 


26837. 


HMX 


.0665 


10230. 


29934. 


EETN 


. 06 25 


9300. 


27231 . 


TETRYL 


. 0585 


8200. 


25755. 


COUP B 


.0607 


8800. 


25722. 


OCTOL 


.0650 


9500 . 


28356. 


Table 


(III-2) . Densities 


of common case 


materials 




Case material 


Density (lb/cu. 


iq) 




Steel 


.IBS 






Aluminum 


. 100 






Uranium 


.688 






Titanium 


. 163 






Lead 


.410 
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CTH 




Figure (III-1). Form of warhead as used by the program 
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a. USER INSTRUCTIONS 

If it is desired to abort the operation of this program 
prematurely, two methods are available. When the program is 
waiting for data entry, press { ENTEfij . When the program is 
not waiting for data entry but is processing, type "HX" and 
press ( ENTERt . Both actions will return the terminal to CMS 
mode . 

When the screen becomes full, or "MORE...." appears in 
the status area, clear the screen by pressing | ALT j and 
1 CLEAR | simultaneously. At several points in the program, 
the user will be directed to "CLEAR SCREEN AND INPUT 0". 

This is to provide proper positioning of the output on the 
screen for ease of reading. If any other symbol than 
"OEFM" should appear in the lower left of the screen, 
press | RESET) aad continue. 

1. Turn the terminal on with the red 

2. When the large "NPS" appears after about 30 seconds, 

press ) RES ET t , then press 1 ENTER} . 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nnnnP", where nnnn is 
ycur 4-digit user number, then press } ENTER 1 . 

4. You will now be asked for your password. Type it in 
(the characters will not appear on the screen), then press 

{ ENTEB 1 . 

5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 

does not, type "GLOBAL TXTLIB FORTMOD2 MOD2EEH" , then press 
lENTERj . 

6. Type "LINK TO xxxxP 191 AS 192 RR", where xxxx is 
the 4-digit user number for the project file, then press 

j ENTER} . 

7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 
then press 1 ENTER j . 



0 



switch. 
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8 . 

9 . 

screen 
10 . 
11 . 



Type "ACCESS 192 fl" and press j ENTER ( 



Press 1 ALT 1 and j CLEAR ) simultaneously to clear 



Type "LBOMB" and press j ENTER j . 

Input the following data as it is reguested on the 
screen by typing the data and then pressing 



ENTER 



Ensure 



that the data is input as either decimal or integer as spe- 
cified. The terms in parenthesis are the program variable 
names. Input the following variables as decimal values: 

Explosive density (XDEN) lb/cu.in 

Explosive Gurney constant (GC) ft/sec 

Explosive detonation velocity (VD) ft/sec 

Case material density (CDEN) lb/cu. ft 

Desired icill radius (MILL) feet 

Warhead diameter (DW) inches 

Warhead length-to-diameter ratio (LTD) 

Target length (LT) feet 

Target mean width (WT) feet 

Target vulnerability, P(k/h) (PKH) 

Target altitude (ALT) 

Target skin thickness (TTH) 

Target skin material (MAT) 



feet 
inches 
1 . 0=Al uminum 
2.0=Fiberglass/Kevlar 
3 . 0=Steel 



12. After entering the above data, you will be presented 
with an initial velocity table built around your desired 
kill radius. The initial velocities will be listed as a 
function of fragment mass and impact velocity. Also pro- 
vided will be the ballistic limit velocities for each of the 
fragment masses. Input the following parameters as decimal 
values: 



Fragment mass (IFMLB (1) ) 
Impact velocity (VHIT(1)) 
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grains 
f t/sec 



13. You will now be presented with a shopping list of 
warheads developed for a range of Pd's. The warheads are 
described by the following quantities: 

Warhead weight in pounds (total weight) 

Case thickness in inches 

Core diameter in inches 

Number of fragments from the warhead 

Number of fragments hitting the target 

Edge length of the fragments in inches 

Input the following parameter as a decimal value: 

Desired probability of a hit given a detonation (PDF) 



14. If you desire to rerun the problem, or want to run a 
new problem, answer the questions appropriately when asked 
by the terminal after the execution of the current problem. 

15. To receive the printout of your encounters, answer 
"nc M to rerunning or restarting the problem when asked by 
the terminal and follow the directions presented on the 
screen. 



16 

press 

17 



Upon completion of the program, type "LOGOFF" and 



ENTEBl . 

Turn the terminal off with the red 



I 

O 



switch. 



C. EXAMPLE PROBLEMS 
1 . Example III-A 

It was desired to build a warhead which would kill a 
typical cruise missile flying at 100 feet altitude. The 
warhead was reguired to have a lethal radius of 50 feet with 
a Pd of 0.98. TNT was selected for the explosive and steel 
was chosen for the case material. The diameter of the mis- 
sile was 13.5 inches. 
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Table (III-3) outlines the input parameters. Table 
(III-4) shows the program output. After the initial veloc- 
ity table was displayed at the terminal, the fragment mass 
was chosen to be 130 grains and the impact velocity was 
selected to be 2500 feet per second. After the Pd table was 
displayed, 0.98 was entered as the kill probability. 

2. Example III-B 

A warhead was required which would kill a typical 
manned aircraft at 30000 feet. A lethal radius of 75 feet 
was specified. The warhead was limited in weight to 50 
pounds and needed to have a core diameter of at least 5 
inches. The diameter of the missile was 10.0 inches. HMX 
was chosen as the explosive and depleted uranium as the case 
material. 

Table (III-5) outlines the input parameters. Table 
(III-6) shows the program output. After the initial veloc- 
ity table was displayed at the terminal, the fragment mass 
was chosen to be 210 grains and the impact velocity was 
selected to be 5000 feet per second. After the Pd table was 
displayed, it was determined that the fragments were too 
large and the fragment mass was then reduced to 100 grains. 
The impact velocity was also reduced to 3000 feet per sec- 
ond. When the Pd table was redisplayed, 0.995 was chosen as 
the desired kill probability. 
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TABLE (I 1 1-3 ) . INPUT DATA FOR EXAMPLE I I I- A 



THE FOLLOWING IS A SUMMARY OF 

01) EXPLOSIVE DENSITY 

02) EXPLOSIVE GURNEY CONSTANT 

03) EXPLOSIVE DETONATION VELOCITY 

04) CASE MATERIAL DENSITY 

05) DESIRED KILL RADIUS 

06) WARHEAD DIAMETE 

07) WARHEAD L ENGTH-TD -DI AMETER RATIO 

08) TARGET LENGTH 

09) TARGET WIDTH 

10) TARGET VULNERABILITY, P(K/H) 

11) TARGET ALTITUDE 

12) TARGET SKIN THICKNESS 

13) TARGET SKIN MATERIAL 



THE INPUT DATA: 



0.05740 


LB/CU. IN 


7600.00 


FT/SEC 


21785.00 


FT/SEC 


0.2830 


LB/CU. IN 


50.0 


FEET 


13.50 


INCHES 


2.50 




36.00 


FEET 


3 .50 


FEET 


0. 150 




100 . 


FEET 


0.060 


INCHES 


STEEL 
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TABLE (I 1 1-5 ) 



INPUT DATA FOR EXAMPLE 1 1 I— B 



THE FOLLOWING IS A SUMMARY OF THE 
01) EXPLOSIVE DENSITY 
021 EXPLOSIVE GURNEY CONSTANT 

03) EXPLOSIVE DETONATION VELOCITY 

04) CASE MATERIAL DENSITY 

05) DESIRED KILL RADIUS 

06) WARHEAD DIAMETE 

07) WARHEAD LENGTH-TO-DI AMETER RATIO 

08) TARGET LENGTH 

09) TARGET- WIDTH 

10) TARGET VULNERABILITY, P(K/H) 

11) TARGET ALTITUDE 

12) TARGET SKIN THICKNESS 

13) TARGET SKIN MATERIAL 



INPUT DATA: 
0. 06650 

10230 .00 

29934.00 
0.6880 

75. 0 
10.00 
2.00 
60.00 
20.00 
0.100 
30000. 
0. 180 
ALUMINUM 



LB/CU.IN 

FT/SEC 

FT/SEC 

LB/CU.IN 

FEET 

INCHES 

FEET 

FEET 

FEET 

INCHES 
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OUTPUT DATA FOR EXAMPLE III- 
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D. PROCEDURAL FLOWCHART 



Input data 

"{ 




Calculate characteristics of 
atmosphere at altitude 


t 




Calculate initial velocities 
for all combinations of preset 
Impact 

velocities and fragment 
weights at the kill radius 




1 


Print table of initial veloci- 
| ties 


l 




Do you want to change the kilx\*_£ , 

^ radius? 

{yes 


| Input new kill radius 








Input desir i fragment mass 
and impact velocity 


<, 



Calculate the critical frag- 
ment spray radius 



the critical radius? 



NO 





Calculate average number of 
hits as a function of 1/R 






Calculate average number of 
hits as a function of 1/R 




, -i , 

[Calculate charge- to-mass ratio] 
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i 

Calculate individual warhead 
solutions for preset values of 
Pd 

Print table of warhead solu- 
tions 

If Pd was restricted due to 
max available warhead radius, 
print 

warning message 




Input desired Pd 



Calculate refined value of 
charge-to-ma ss ratio 

Determine warhead solution for 
specified Pd 




Print warhead solution 
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IV 



PROPULSION MOTOR SIZING 



A. DESCRIPTION AND ORIGIN OF PROGRAM 

This program provides a method for the preliminary siz- 
ing of a solid propellant rocket motor for a boost-sustain 
trajectory of a tactical missile. The analytical method was 
developed by Redmon (Ref. 1] and was expanded with the addi- 
tion of material from Platzek [Ref. 2] and Hill (Ref. 3]. 

The program was written for this thesis in FORTRAN IV for 
use on the IBM 370 computer system. Essentially, the analy- 
sis consists of sizing the booster motor from differential 
velocity and acceleration requirements with limitations 
imposed by the physical dimensions of the missile. The 
booster is at all times considered to be a core-burning 
motor. The sustainer motor calculations are controlled by 
the maximum range sp-ecified by the user and by the solution 
of the booster motor. The sustainer can be either a core- 
burning or an end-burning motor. 

The rocket motor configuration is assumed to be either 
an integral booster- sustainer motor as shown in Figure 
(I V— 1 ) or a staged motor as in Figure (IV-2) . The booster 
and the sustainer always burn exclusively, or, in other 
words, one is not burning while the other one is. The 
nozzle half angle is specified by the user consistant with 
space available in the missile. If a staged motor is uti- 
lized, both nozzles will have the same half angle. 

The booster calculations start by determining the amount 
of thrust needed to boost the total weight of the missile to 
its cruise velocity at the prescribed acceleration. From 
this, the necessary amount of propellant is estimated and 
the process is iterated to account for the decreasing mass 
situation. The chamber pressure is initially estimated by 
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Figurre (IV-1). Integral booster-snstainer notnr 
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PROPELLANT 




Figure (IV-2) . Staged motor 
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BOOSTER 1 BOOSTER * SUSTAINER * SUSTAINER 

NOZZLE NOZZLE 



minimizing the motor weight to propellant specific impulse 
ratio, as presented by Eedmon [Ref. 1], The next step of 
the program is to size the nozzle using the now known ini- 
tial pressure ratio. Once the ideal nozzle is developed, 
the chamber pressure is raised or lowered as necessary in an 
attempt to cause the nozzle exit diameter to exactly match 
the booster diameter. However, if the integral motor option 
is used, the program will drive the nozzle as small as it 
can without violating one of the following limits in order 
to increase the probability that the sustainer will operate 
properly. The iterative process has two limits: an abso- 

lute maximum of 2000 PSI chamber pressure [Ref. 2] and a 
minimum of 1000 PSI if the pressure had previously been 
higher. The solution can be less than 1000 PSI if the pres- 
sure remained below that level throughout the problem. Also 
limited is the exit pressure to ambient pressure ratio. At 
the low end, it is limited to 0.4 to prevent flow separation 
in the nozzle. At the high end, it is limited to 2.5 to 
prevent excessive underexpansion and loss of physical real- 
ity in the program results [Ref. 4], The burn rate is 
initialized at a potential maximum of 1.25 inches/second 
(Ref. 2} and is allowed to decrease to arrive at a compati- 
ble burn area and web thickness combination. 

The sustainer motor, in the integral motor case, is vir- 
tually a continuation of the booster solution. The initial 
thrust requirement is determined by increasing the cruise 
speed drag to account for speed loss through maneuvers. It 
is also then refined for weight change if any climbing or 
diving is required to reach the target altitude. The nozzle 
is the same one as developed for the booster except that the 
throat area will be expanded as a result of the erosion 
effect. From the new area ratio, a pressure ratio is deter- 
mined. The chamber pressure and thrust coefficient are then 
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iterated until a steady chamber pressure evolves to provide 
the required thrust. If at any time it drops below 125 ESI, 
the program will stop since this is considered a minimum 
chamber pressure for proper propellant combustion. The exit 
pressure to ambient pressure ratio remains subject to the 
same restrictions. The burn rate starts at 0.45 inches/sec- 
ond £Ref. 2] and is decreased to provide an acceptable web 
thickness and burn area. The solution can be either an end 
burner or a core burner, depending on the burn area 
required. 

The sustainer for the staged motor is solved in essen- 
tially the same manner as the booster. The two major excep- 
tions are that it can produce an end burner if the burn area 
is small enough and the thrust required is based on the 
cruise drag and not the velocity to be gained. The chamber 
pressure is limited to an absolute maximum of 800 PSI and a 
minimum of 250 PSI if the pressure had previously been 
higher [Her. 21. The same pressure ratio restrictions apply 
and the burn rate starts at 0. 45 inches/second. 

Other motor-nozzle combinations can be created from the 
results of this program on a first-order approximation 
basis. Figures (IV- 3) and (IV-4) illustrate two possible 
methods for separating the motor nozzles without resorting 
to staging. Although these two methods will probably have 
the same weight disadvantage that staging does, they can be 
packaged in to a smaller volume of space. A hint for "crea- 
tive construction" is to rerun the problem after sizing the 
motors using a smaller missile diameter to force the nozzles 
to a smaller size. 
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Figure (IV- 3) . Separate nozzles (nonstaged moto 
The values are obtained from the staged motor opt 
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Figure (IV-4) . Concentric nozzles (nonstaged actors) 
The values are obtained from the staged motor option. 
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B. USES INSTR UCTIONS 

If it is desired to abort the operation of this program 
prematurely, two methods are available. When the program is 
waiting for data entry, press (enter! . when the program is 
not waiting for data entry but is processing, type "HX" and 
press | ENTER! ♦ Both actions will return the terminal to CMS 
mode . 

When the screen becomes full, or "MORE...." appears in 
the status area, clear the screen by pressing | ALT I and 
I CLEAR | simultaneously. At several points in the program, 
the user will be directed to "CLEAR SCREEN AND INPUT 0". 

This is to provide proper positioning of the output on the 
screen for ease of reading. If any other symbol than 

" should appear in the lower left of the screen, 
press 1 RES El | and continue. 

1. Turn the terminal on with the red 

2. When the large " NPS" appears after about 30 seconds, 

press | BESET 1 , then press | ENTER| . 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nnnnP", where nnnn is 
your 4-digit user number, then press | ENTER | . 

4. You will now be asked for your password. Type it in 
(the characters will not appear on the screen), then press 

1 ENTER f . 

5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 

does not, type "GLOBAL TXTLIB F0RTM0D2 MOD2EEH" , then press 

1 ENTER) . 

6. Type "LINK TO xxxxP 191 AS 192 R B", where xxxx is 
the 4-digit user number for the project file, then press 

1 ENTER! . 

7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen), 
then press j EN TER | . 



switch. 
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8. Type "ACCESS 192 B" and press j ENTER | . 

9. Press ( ALT | and | CLEAR j simultaneously to clear 



screen. 

10. Type "LPROP" and press j EN TER{ . 

11. Input the following data as it is requested on the 
screen by typing the data and then pressing 



ENTER 



Ensure 



that the data is input as either decimal or integer as spe- 
cified. The terms in parenthesis are the program variable 
names. 

The following variables are required inputs for both 
motor option problems. 



Motor option 


(IMOTO R) 


0=inregral motors 
1=staged motors 


Launch altitude 


(LALT) 


feet 


Launch weight 


(ML) 


pounds 


Launch velocity 


(VBI) 


feet/second 


Launch elevation angle 


(ELB) 


degrees 


Boost acceleration 


(a) 


gravities 


Cruise velocity 


(V3F) 


feet/second 


Cruise velocity drag 


(DRAGS) 


pounds 


Maximum range to target 


(R) 


nautical miles 


Maximum target altitude 


(TALT) 


feet 


Booster propellant 
specific impulse 


(ISP3) 


seconds 


Booster propellant 
density 


(DENSB) 


pounds/cu. inch 


Booster exhaust 
specific heat ratio 


(GB) 




Sustainer propellant 
specific impulse 


(ISPS) 


seconds 


Sustainer propellant 
density 


(DENSS) 


pounds/cu. inch 


Sustainer exhaust 
specific heat ratio 


(GS) 




Nozzle half angle 


(ALN) 


degrees 
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The following variables are required inputs for the 
integral motors option only. 



Nozzle design altitude (ALTBN) 



Nozzle erosion rate (EH) 

Missile diameter (DB) 

Case yield strength (YIELD) 

Case density (DENSC) 



feet 

inches/second 

inches 

PSI 

pounds/cu. inch 



The following variables 
staged motors option only. 

Booster design altitude 

Booster diameter 

Booster case 
yield strength 

Booster case density 

Sustain er design 
altitude 

Sustainer diameter 

Sustainer case 
yield strength 

Sustainer case density 



are reguired inputs for the 



(ALTBN) 


feet 


(DB) 


inches 


(YIELD) 


PSI 


(DENSC) 


pounds/cu. inch 


( ALTS N) 


feet 


(D) 


inches 


(YIELD) 


PSI 


(DENSCS) 


pounds/cu. inch 



12. This program will cue the user when the input param- 
eters have dictated a scenario which either cannot be 
achieved within reality or produce less than optimum 
requirements on the propulsion system of the missile. They 
are not definitive and are only intended to make the user 
aware of a situation which may need correction. The follow- 
ing is a list of available cue messages with short 
definitions. 

"SUSTAINER NOT CALCULATED SINCE THE BOOSTER BURNOUT 
RANGE EXCEEDS THE DESIGN RANGE." This can result from 
entering an extremely short range for the missile, or it 
can be caused by a very low acceleration requirement. 

"SUSTAINER NOT CALCULATED: THE BOOSTER NOZZLE DESIGN 

PREVENTS SUSTAINER OPERATION. RECOMMEND STAGING OR 
INDEPENDENT NOZZLES." This occurs only when using the 
integral motors option. Then scenario described to the 
program can cause the booster nozzle to be too large to 
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maintain the proper chamber pressures when the motor has 
shifted to sustamer operation. This usually occurs 
when a large acceleration is demanded but the thrust 
required fcr the cruise trajectory is small. 

"BOOSTER BURN RATE MAS ADJUSTED TO CORRESPOND WITH THE 
W£E THICKNESS." This occurs quite often and simply 
indicates that the burn rate was decreased from its 
potential physical maximum of 1.25 inches/second. 

"THE BOOSTER NOZZLE DESIGN WAS RESTRICTED DUE TO THE 
MISSILE DIAMETER." and "THE S USTAINER NOZZLE DESIGN WAS 
RESTRICTED DUE TO THE MISSILE DIAMETER." The nozzle was 
not able to be designed for optimum pressure ratios at 
the mid point of the boost trajectory. Usually, the 
exit diameter is solved larger than the missile diameter 
and is subsequently reduced to fit. 

"BOOSTER NOZZLE DESIGN IS NOT OPTIMUM DUE TO EXCESSIVE 
BOOSTER CHAMBER PRESSURES," and " SUSTAINER NOZZLE DESIGN 
IS NOT OPTIMUM DUE TO EXCESSIVE SUSTAINER CHAMBER 
PRESSURES." If the nozzle cannot be downsized without 
exceeding pressure thresholds (2000 PSI for the booster 
and 800 PSI for the sustainer) . the chamber pressure is 
held just below the pressure threshold and the nozzle 
area ratio will be adjusted to allow the nozzle to fit 
in the missile. 

"THE SUSTAINER MOTOR HAS A CORE-BURNING GRAIN." The 
required burn area for the sustainer was too large too 

E ermit an end burning grain with a properly realistic 
urn rate. 

"SUSTAINER BURN RATE WAS ADJUSTED TO CORRESPOND WITH THE 
WEB THICKNESS." This indicates the sustainer burn rate 
was lowered from a potential maximum of 0.45 inches/sec- 
ond to provide a proper web thickness. 

"THE SUSTAINER MOTOR HAS AN END BURNING GRAIN." The 
required burn area for the sustainer was small enough to 
permit an end burning grain. The burn rate is then 
adjusted to correspond with the nonreduceable burn area. 

"REESTIMATION OF LAUNCH WEIGHT IS REQUIRED FQR THESE 
MISSILE PERFORMANCE PARAMETERS." The scenario described 
to the program produced a motor whose weight is either 
larger than 75ft of rhe total or less than 25ft of the 
total. 

"ENLARGEMENT OF DIAMETER IS RECOMMENDED DUE TO A VERY 
HIGH LENGT H-TO- DIAMETER RATIO FOR THE MOTOR." This cue 
indicates the length-to-diamet er ratio is greater than 
15* Other components of the missile will make the mis- 
sile's overall length-to-diameter ratio even larger. 



13. Immediately after completion of the solution, the 
pregram will ask if you want to receive a hardcopy printout 
of that particular solution. A "yes" answer stores that 
solution in a file for retrieval by the user when he fin- 
ishes with the program. 
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14. If you desire to rerun the problem, or want to run a 
new problem, answer the questions appropriately when asked 
by the terminal after the execution of the current problem. 

15. To receive the printout of your encounters, answer 
"no" to rerunning or restarting the problem when asked by 
the terminal and follow rhe directions presented on the 
screen. 



16. Upon completion of the program, type "LOGOFF" and 



press | ENTEB( . 

17. Turn the terminal off with the red 



I 

O 



switch. 



C. EXAMPLE PROBLEMS 

1 . Example IV-A . Integral motors, common nozzle 

The following parameters are input for the integral 
motor example: 



Launch altitude 


35.0 


feet 


Launch weight 


1000.0 


pounds 


Launch velocity 


0.0 


feet/second 


Launch angle 


60.0 


degrees 


Average acceleration 


30.0 


g *s 


Cruise velocity 


4000. 0 


feet/second 


Drag at cruise velocity 


1500.0 


pounds 


Maximum range 


20.0 


miles 


Final (target) altitude 


50000 


feet 


Booster propellant 
specific impulse 


260.0 


seconds 


Booster propellant 
density 


0.075 


lbs/cu. inch 


Booster exhaust specific 
heat ratio 


1. 244 




Sustainer propellant 
specific impulse 


210.0 


seconds 


Sustainer propellant 
density 


0. 065 


lbs/cu. inch 


Sustainer exhaust 
specific heat ratio 


1.270 




Nozzle half angle 


20.0 


degrees 


Nozzle design altitude 


0.0 


feet 
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TABLE <IV-1K OUTPUT OF EXAMPLE IV- 
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Nozzle erosion rate 
Missile diameter 



0.001 inches/second 
10.0 inches 



yield strength of 

case material 180000.0 PSI 

Density of case 

material 0.2662 lbs/cu.inch 

The solution for this problem is presented in Table 

(XV- 1) . 

2. Sxample IV-B . Staged motors, separate nozzles 



following are input for 


the staged motor prblem 


Launch altitude 


35.0 


feet 


Launch weight 


2200.0 


pounds 


Launch velocity 


0.0 


feet/second 


Launch angle 


30.0 


degrees 


Average acceleration 


25.0 


g ' s 


Cruise velocity 


2200.0 


feet/second 


Drag at cruise velocity 


1000.0 


pounds 


Maximum range 


50.0 


miles 


Pinal (target) altitude 


75000 


feet 


Booster propellant 
specific impulse 


250.0 


seconds 


Booster propellant 
density 


.0647 


lbs/cu.inch 


Booster exhaust 
specific heat ratio 


1.225 




Sustainer propellant 
specific impulse 


205.0 


seconds 


Sustainer propellant 
density 


. 0625 


lbs/cu. inch 


Sustainer exhaust 
specific heat ratio 


1.257 




Nozzle half angle 


15.0 


degrees 


Booster design altitude 


0.0 


feet 


Booster diameter 


14.5 


inches 


lield strength of 
booster case 


180000.0 


PSI 


Density of booster 
case material 


0.2662 


lbs/cu. inch 


Sustainer design altitude 0.0 


feet 
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Sustainer diameter 



14.5 inches 



Yield strength of 

sustainer case 180000.0 PSI 

Density of sustainer 

case material 0.2662 lbs/cu.inch 

The output for this example is provided in Table 



(IV-2) . 
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PROCEDURAL FLOWCHART 



Inpat data 



I 



Calculate average drag of 
boost phase based on cruise 
speed drag and initial and 
final speeds 



I 



Estimate thrust required for 
boost and calculate booster 
burn time 



I 



Determine averaae altitude 
during boost and calculate 
ambient pressure 



Refine thrust required for 
boost and determine weight of 
booster propellant 



Determine mass flow rata 



Calculate pressure ratio based 
on optimum hamber pressure 
to case weight relationship 



PLIM=2000 PSI 



Calculate thrust coefficient, 
throat diameter and exit 
dia meter 



<■ 



Is exit diameter smaller tha 
base diameter? 



I 



r\ 



MO 



Increase chamber pressure 

"1 
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i 
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{ 

Set burn rate equal to 1.25 
and calculate resulting web 
thickness 



- J 


Determine burn area and 
| required web thickness 


1 

< 

1 


- 1 - 

Are the two web thickness 
values equal? 

r--~ " 


Lower burn rate 


I 



Calculate di nensions of the 
booster and nozzle 



Is this a staged motor? y 



wo 



Calculate reguired sustainer 
burn time 



Determine average altitude 
during cruise ana calculate 
ambient pressure 



Calculate thrust required to 
offset cruise drag and 
calculate weight of sustainer 
pru oellant 



I 



Determine mass flow rate and 
nozzle area ratio 



:::::::::: it::::::::::::: 

[Calculate Mach of the exhaust 

"t 
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Calculate nozzle pressure 
ratio and estimated thrust 
coefficient 



_»{ Calculate sustainer chamber 
pressure 



<! 



Is chamber pressure less tha 



^25 PSI? 



3 






HO 



Refine thrust coefficient 



YES 



estimated thrust coefficient? 

" 1 


1 

L 


j- Set them equal 




j Set burn rate equal to 0.45 j 


r 



J 


[ Determine required burn area | 


n 

< 


1 ^ 


Is the burn area less than the 
motor cross sectional area? 




Calculate web thickness | 




i 


s' Is the web thickness less than 


86% of the motor radius? 




1 MO 


1 


| Lower the burn rate 






Yes 




Set burn ar a equal to motor 
cross-section area and lower 
burn rate to correspond 



I 



Establish dimensions of 
tainer motor 



sus- 
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(Calculate sustainer burn time 



Determine average altitude 
during cruise and calculate 
ambient pressure 



Calculate thrust required to 
offset cruise drag and 
calculate weight of sustainer 
propellant 



Determine mass flow rate 

Calculate pressure ratio based 
on optimum hamber pressure 
to case weight relationship 

PLIM=800 PSI 

Calculate thrust coefficient, 
throat dianeter and exit 
diameter 



Mo 



< Is exit diameter smaller than 
base diameter? 

Increase chamber oressure 



yes 



< 



Is chamber Dressure larger 
than PLIM? 



J_Yes 



> 



Lower chamber pressure 
Increase exit pressure 



I 



< Is (exit pressure/ambient 
pressure) greater than 0.4? 
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7. AERODYNAMIC COEFFICIENTS 



A. DESCRIPTION AND ORIGIN 

This program is the current edition of a program which 
originated at the Naval Ship Research and Development Center 
in 1971 [Ref. 5J. The program was written in FORTRAN II for 
use on the IBM 7090 digital computer. It was developed as a 
method for predicting the static, longitudenal aerodynamic 
characteristics of typical missile configurations with the 
ccntrcl surfaces in a •'plus' 1 attitude. The original program 
computed the aerodynamic characteristics for missiles oper- 
ating at angles of attack up to 180 degrees. The effects of 
control surface deflections for all modes of aerodynamic 
ccntrcl are taken into account. The method was based on the 
then well known linear, nonlinear crossflow and slender body 
theories with empirical modifications to provide the high 
angle of attack capability. 

The program was modified and, presumably, updated in 
1 S74 by F. A. Kuster, Jr., of the Naval Air Development Cen- 
ter. In 1980, the program was modified for use on the Naval 
Postgraduate School IBM 360 computer system by D. Redmon 
[Ref. II. The current version of the program was modified 
fcr use on the new Naval Postgraduate School IBM 370 compu- 
ter system. It has been expanded to provide graphical pre- 
sentation of the output data. 

It must be emphasized at this point that the current 
program edition is not in a completely finished state. 
Somewhere in the history of the program after its initial 
establishment on the IBM 360, errors were introduced during 
the modifications. At present, these errors do not prevent 
the use of the program and the output data is considered to 
be correct for trend-observance purposes. 
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Specifically, the program does not produce any drag-rise 
phenomena for either the wings or the tails when Cd is 
observed as a function of Mach number. Additionally, the 
decline of the drag coefficient above Mach 1.0 is not smooth 
or as prolonged as is found experimentally. It is very pro- 
bable that these two failings of the program are linked to a 
common error inserted accidentally in the process of tailor- 
ing the program for use on the IBM 370. In order to 
temporarily smooth over the graphical discontinuities, expo- 
nential decay functions were inserted. They are clearly 
marked in the the program listing for removal when the pro- 
gram is corrected. 

The input to the program is composed of a detailed list- 
ing of the dimensions of the missile to be analyzed. The 
current version of the program will consider a missile which 
has four symmetrical wings and four symmetrical tails. The 
missile may be either canard or tail equipped and either 
wing or canard or tail controlled. The program assumes the 
control surface is the tail, however, the input data is 
"mislabelled" to produce the proper configuration. For ins- 
tance, if the missile is a wing control missile, the wing 
data is input as the tail and the tail data as the wing. 

For a canard controlled missile, the canard data is input as 
the tail. Figures (7-1) and ( V— 3) show two typical missile 
configurations and where the input parameters for the pro- 
gram are obtained. 

B. USEE INSTRUCTIONS 

If it is desired to abort the operation of this program 
prematurely, two methods are available. When the program is 
waiting for data entry, press | BNTBR| . When the program is 
not waiting for data entry but is processing, type "HX" and 
press j ENTER | . Both actions will return the terminal to CMS 
mode. 
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When the screen becomes full, or " MORE. . . . " appears in 
the status area, clear the screen by pressing j ALTj and 
J CLEAR t simultaneously. At several points in the program, 
the user will be directed to •• CLEAR SCREEN AND INPUT O'*. 

This is to provide proper positioning of the output on the 
screen for ease of reading. If any other symbol than 

should appear in the lower left of the screen, 
press j RESETj and continue. 

1. Turn the terminal on with the red 

2. When the large "NPS" appears after about 30 seconds, 

press 1 RESET} , then press \ ENT ER| . 

3. When "CP READ" appears in the status area on the 
lower right of the screen, type "L nnnnP", where nnnn is 
ycur 4-digit user number, then press ( ENTER j . 

4. You will now be asked for your password. Type it in 
(the characters will not appear on the screen), then press 

| ENTER j • 

5. Your personal file must contain a PROFILE EXEC rout- 
ine with the appropriate Fortran GLOBAL statement. If it 

does not, type "GLOBAL TXTLI3 FORTMOD2 MOD2EEH" , then press 
jENTEBf . 

6. Type "LINK TO xxxxP 191 AS 192 RR", where xxxx is 
the 4-digit user number for the project file, then press 

[ ENTER J . 

7. You will now be asked for the project file password. 
Type it in (the characters will not appear on the screen). 



then press | ENTER | . 



8. 


Type " 


ACCESS 192 B" and press | ENTER | . 


9. 


Press 


) ALT) and j CLEAR] simultaneously to clear 


screen. 






10. 


Type " 


LAEROI " and press IenterI. 


11 . 


Input 


the following data as it is reguested on the 



screen by typing the data and then pressing | ENTER} . Ensure 



switch. 
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that the data is input as either decimal or integer as spe 
cified. The terms in parenthesis below are the program 
variable names. 



Input the following as integer values unless otherwise 
noted. The integers must be two digit integers (three=03) . 



Control constant (ICSC) 



01- Tail control 

02- Wing control 

03- Canard control 



Hose shape (INOSE) 



01- Ellipsoid 

02- 0give 

03- Cone 



Humber of control deflections (IDT) Less than 11 

You will now be asked for the control deflections in 
degrees, decimal values. 



Humber of Hach numbers (IM) Less than 25 

You will now be asked for the Hach numbers, decimal 
values. Each Mach number will produce a separate 
table and plot of output data. 

Humber of angles of attack (IAL) Less than 25 

You will now be asked for the angles of attack in 
degrees, decimal values. 



Humber of configurations { HBOD Y) Ho restrictions 

Each configuration will restart the program. Only 
the last configuration will produce the written 
output. 

Wing planform (ISWPW) 01-Not delta 

02-Delta 



Wing position (IAFBW) 

Wing sweep constant (ISWEPW) 

Humber of wings (NWIHG) 

Tail planform (ISWPT) 



00- Body after wing 

01- No body after wing 

00- Delta planform or 
Unswept leading 
edge 

01- Swept leading edge 
04 

01- Not delta 

02- Delta 



Tail position (IAFBT) 

Tail sweep constant (ISWEPT) 

Humber of tails (NTAIL) 



00- No body after tail 

01- Body after tail 

00- Delta planform or 
Unswept leading 
edge 

01- Swept leading edge 
04 - 



Input the following values as decimal numbers: 



Wing tip-to-chord ratio (XLAMW ) 



Wing leading edge sweep (CLAMW) Degrees 
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Wing span including body (BW) 

Wing root chord (CROOTW) 

Exposed wing area, 2 panels ( S W) 

Wing mean geometric chord (XMACW) 

Distance from nose to 

wing leading edge (XWING) 

Wing thickness-to-chord ratio (TOVCW) 

Tail tip-to-chord ratio (XLAHT) 

Tail leading edge sweep (CLAMT) 

Tail span including body (3T) 

Tail root chord (CROOTT) 

Exposed tail area, 2 panels (ST) 

Tail mean geometric chord (XMACT) 

Distance from nose to 

tail leading edge (XTAIL) 



Tail thickness-to-chord ratio (TOVCT) 
Altitude (HT) 


Feet 


Body diameter (D) 


Fee t 


Missile length (XL) 


Feet 


Nose length (XLNOSE) 


Feet 


Distance from nose to CG (XCG) 


Feet 


Reference area (AREA) 


Square feet 


Reference length (XREF) 


Feet 


Engine code (ENGINE) 


0. 0-Turbofan 
1 . 0-Rocket 


Inlet code (ENLET) 


0 . 0-Flush 
1 . O-Extended 


Boat tail angle (BETA) 


Degrees 


Base diameter (DBASE) 


Feet 


Nozzle exit diameter (DJET) 


Feet 


Boat tail length (XLABOD) 


Feet 


Protuberance drag (CDPROT) 


(Coefficient value) 


If comparing with experimental 
values, Reynolds number (REFT) 


(Dimensionless) 



Degrees 

Feet 

Feet 

Square feet 
Feet 

Feet 



Feet 

Feet 

Square feet 
Feet 

Feet 
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12. If you desire to rerun the problem, or want to run a 
new problem, answer the questions appropriately when asked 
by the terminal after the execution of the current problem. 

13. To receive the printout and plot of your encounters, 
answer "no" to rerunning or restarting the problem when 
asked by the terminal and follow the directions presented on 
the screen. 

14. Upon completion of the program, type "LOGOFF" and 



press J ENTER} . 

15. Turn the terminal off 



with the red 




switch. 



C. EXAMPLE PEOBLEMS 

Table (V-5) identifies the output variables as they 
appear in the output tables. 

1 . Example V-A . Tail control missile 

Figure (V-1) illustrates the missile used in this 
example. The dimensions for this missile and ether input 
parameters are contained in Table (V-1) . The output is 
shown in Table (V-2) and Figure (V-2) . 

2. Example V-II . Canard control missile 

Figure (V-3) illustrates the conard configuration 
missile used in this example. The input data is contained 
in Table (V-3) . The output is displayed in Table (V-4) and 
Figure (V-4) . 
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Figure (V- 1) . 



Tail control missile as used in Example V-A 
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T HE FCLLCWING T A 0L E CONTAINS THE INPUT DATA FOR 
EXAMPLE V-A. TAIL CONTROL MISSILE 



1) (ICSC) CCNTRCL CONSTANT; I = TA I L ♦ 2=WING, 3 = C AN ARD: 

2) ( INCSE) NCSE jH AP E ; 1=ELLIPSE, 2-OGIVE, 3 = CGNE ; 

3) (IDT) NUM.3ER OF CLN T R l L DEFLECTIONS: 

4) (IP) NUMBER CF 4ACH NUMBERS: 

5) ( I A L ) NUMBER CF ANGLES OF ATTACK: 

6) (N8GOY) NUMBER OF CONFIGURATIONS: 

7) (ISwPW) L -NUN-DEL T A W i NG , 2 -DEL T A WING: 

8) (IAFBW) 3-Nu BODY AFTER WING, 1-BUOY AFTER WING: 

9) (IS*EPW) WING SrtEtP CONSTANT i I F DELTA-0) 

0=UN SWEPT LEADING EDGE, I=SWEPT LEADING 

10) (NWING) NUMBER GF WINGS: 

11) ( I S rt P T ) 1=NQN-DELTA TAIL, 2 =DE L TA TAIL: 

12) (IAFBT) 0- NG BODY AFTER TAIL, 1-BGDY AFTER TAIL: 

13) (ISwEPT) TAIL S w E r P CONSTANT (IF CELT A=0 ) 

0= UNSW EPT LEADING EDGE, l^SWE^T LEADI NG 

14) (NTAIL) NUMBER OF TAILS: 

L 5 ) (XLAMw) T IP- lu-CHCR.) RATIO OF wING: 

16) (CLAMW) WlNo LEAu I No EDGE SWEEP ( DEGREES) : 

17) (Brt) WING SPA N f INCLUDING BODY: 

10) (CRGUTW) aING ROOT CHORD (AT BODY JJNCTION): 

19) (Srv) EXPOSED WING AREA (T,;Q PANELS): 

20 ( X M A C w ) WlNo MEAN GtOETRIC CHORD: 

21) (XWING) Cl STANCE FRuM NOSE TO WING LE: 

22) (TCVCW) WIN, TH I C.\Nc S $ TO CHORD RATIO: 

23) (XLAMT) TIP-1G-CHLRJ PATIO CF TAIL: 

24) (CLAMT) TAIL LEADING cOGE SWEEP (DEGREES): 

25) (OT) TAIL SPAN, l.NCLUCING BODY: 

26) (CRCOTT) TAIL ROOT CHORD : 

27) (ST) EXPOSED TAIL AREA ( T*0 PANELS): 

(XMACT) TAIL MEAN GEOMETRIC CHORD: 

( XT A I L ) DISTANCE F R J M NOSE TO TAIL LE: 

( 10 VC T ) TAIL IHICaNESS TO CHORD RATIO: 

31) (HT) ALTITUDE: 

32) (C) BODY DIAMETER: 

33) (XL) MISSILE LENGTH: 

34) ( XLNGSE ) NO LENGTH: 

35) ( X C G ) DISTANCE TO CG LOCATION FROM NOSE: 

36) (AREA) REFERENCE AREA: 

37) (XREF) REFERENCE LENGTH: 

38) (ENGINE) ENGINE; C • )= TURBO FAN , i.O=POCKET: 

39) (ENLET) INLcT; U.O=rLoSH, L •0=EXTENDED: 

4G ) (BETA) OCAT-TAlL ANGLc (DEGREES): 

41) (DBASE) BASE DIAMETER: 

42) (DJGT) NOZZLE cKIJ DIAMETER: 

43) (XLA6CC) duAT-TAlL LENGTH: 

44) (CDPRCT) PROTUBERANCE ORAG: 



28) 

29 ) 

30) 



EDGE: 



EDGE: 



1 

3 
5 

1 

11 

1 

2 

1 

0 

4 
L 
1 



1 
4 

0.0 

77.000 
1.800 
2.960 
2. 072 
1.973 
1.902 
0.030 
0 .609 
15. OOC 
1.700 
0.400 
0.418 
0.328 
5.42C 
0.076 
30D00 .000 
0.400 
6.000 
0.750 
3.200 
0.127 
0.400 
1.0 
0.0 
0.0 
0.400 
0.250 
0.0 
0.0 



MACH 


2. SCO 
0.0 


0.0 

c.o 


oo 

• . 

oo 


0.0 

0.0 


0.0 
D .0 


0.0 

0.0 


0. 0 
0.0 


0. 0 
0.0 


0 . c 
0.0 


oo 

• « 

oo 


o.c 

o.c 


OELTA 


o.c 

o.c 


4.00 

c.o 


3.00 

0.0 


12.00 

0.0 


16.00 

0.0 


0.0 

0.0 


0.0 

0.0 


0. 0 
0. 0 


o.c 

o.c 


0.0 
0. 0 


0.0 

0.0 


ALPHA 


o.c 

o.c 


L.OO 

0.0 


2.00 

0.0 


3.00 

0.0 


4.00 

0.0 


5.00 

0.0 


6.00 

0.0 


7. 00 
0.0 


8.00 

0.0 


9. 00 
0.0 


10 . oc 
0. c 



Table (7-1). Inpat data for Example V-A 
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Table (V-2) . Dutput of Example V-A 
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Figure ( V— 2) . Output data plot for Example V- A 
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Figure (V-3) . 



Canard control missile for Example V-B 



THE FCLLCVsING TABLE CONTAINS THE INPUT DATA FOR 
EXAMPLE V- 6 . CANARD CONTROL MISSILE 



1 ) 

2 ) 

3 ) 

4 ) 

5 ) 

6 ) 

7 ) 

8 ) 
9 ) 

10 ). 

ID 

12 ) 

13 ) 

14 ) 
13 ) 
16 ) 

17 ) 

18 ) 

19 ) 

20 ) 
21 ) 
22 ) 

23 ) 

24 ) 

25 ) 

26 ) 

27 ) 

28 ) 
29 ) 

Jfl 

32 ) 

33 ) 

34 ) 

35 ) 

36 ) 

37 ) 
3 d) 

39 ) 

40 ) 

41 ) 

42 ) 

43 ) 

44 ) 



( I C SC ) control constant; 1 = TAIL, 2 =WING, 3 =C ana rd 

(INGSE) NO Sr SHAPE; UELLI^SE, 2 =UG I V E » 3 = CJNE; 
(IDT) N JM 3 E< UF CONTROL DEF LEC T I QMS : 

(IN) NUMBER CF MACH NUMBERS: 

( I A L ) NUMBER CF ANGLES UF ATTACK: 

(NBLDY) NUMBER JF CONFIGURATIONS: 

(ISWPW) 1 = .NON-DEL T A WING, 2 =Dc L TA WING: 

(IAFBW) C = NO -sUJY AFTER WING, L -BOD Y AFTER * ING : 
(ISGEPW) MNG SWEEP CONSTANT (IE DEL T A= 0 ) 

0 - L N S N E PT LEADING EDGE, 1 =SWEPT 
(NWINC-) NUMBER UF MNbS: 

(ISaPT) 1 *NGN- 0 ELTA TAIL, 2 = DEL T A TAIL: 

(IAF 3 T) G = Nw BCD Y AFTER T A I L , 1 = 8 GD Y AFTER 
(ISwEPT) tail swEt-P Constant (if jelta= 0 ) 

0 = UN S WE ^ T LEADING EOGE, l=SWEPT 
( NT A I L ) NUMBER OF TAILS: 

(XLAMW) TIP-To-CHbRJ RATIO CF WING: 

(LLAyfc) ViINj LEADiNi EDwE S WEEP ( DEGRE ES ) : 

( 3 M in IN G SPAN, INCLUDING BUOY: 

(CRCCT.n) *1 JG ROOT CHORD (AT BODY JUNCTION) 



LEADING EDGE: 



TAIL: 

LEADING cDGE : 



(SW) EXPOSE 3 
( XMACW) 

( Xh ING ) 

( TQVCU) 

( XLAMT ) 

(CL AMT) 



WING ARE-* (TrtU PANELS) 

MNG IcAN GEOMETRIC CHORD: 

DISTANCE FROM NOSE TO WiNG LE: 
a ING I H I C K N w S S TO CHORD RATIO: 
TIP- 70 -CH„RJ RATIO OF TAIL: 

TAIL LEADING EDGE SWEEP (DEGREES): 
(BT) TAIL SPAN, INCLUDING BODY: 

(CRCCTT) TAIL ROOT CHORD: 

(ST) EXPOS E) TAIL ARE * (TWO PANELS): 
(XMACT) T ML McAN GEGMETRIu CHORO: 

(XTAIL) DISTANCE FROM NOSE TO TAIL LE: 
(TOVCT) TAIL THICKNESS TO 3 H 0 RD RATIO: 

( HT ) ALTITUDE: 

(0) BCUY DUMB TER: 

(XL) MISSILE LENGTH: 

(XLNCSE) NOSE LENuTT: 

(XCG) OISTAIC*. TU Co LOCATION FROM NOSE: 
(AREA) REFERENCE ARcA: 

(XREF) REFERENCE LENGTH: 

(ENGINE) ENGINE; 0 . J = I URBUF AN , L . 0 = 9 OCKET : 
(ENLET) INLET; 0 . 0 =FLJSH, l . 0 = E X TEND E D : 
(BETA) BGAT-TML AN iLE (DECREES): 

( L B A S E ) B A Sc DIAMETER: 

(CJET) NCZZlE EXIT DIAMETER: 

(XLAOCD) BOAT-TAIL LENGTH: 

(CCPRLT) PROTUBERANCE DRAG: 



3 
2 
5 
1 

11 

1 

1 

1 

1 

4 

l 

i 



1 
4 

0.314 

5 T .500 

3.600 
A. 140 
8.070 

2 . 522 
6.060 
0.0 50 
0 . 541 

32.000 
2.3 00 
0.6 14 
0.568 
0.473 
2.700 
0.085 

70300.000 

3. BOG 

12. 000 

2.20 0 

7.600 
0.503 
0 . 3 00 
1 . 3 

0 . 0 
5.000 
0.600 
0 . 600 
1 . 20 C 
0.0 



MACH 


3 . 5 CO 

o.c 


C.C 0 

0.0 0 


. 0 
. 0 


0.0 

0.0 


0.0 

0.0 


0.0 

0.0 


0. 0 
0 . 0 


0 . 0 
0.0 


0.0 

0.0 


G .0 
0.0 


0.0 

0.0 


DELTA 


c.c 

o.c 


4 . GO 
0.0 


d .00 

3.0 


12.00 

0.0 


16.00 

0.0 


0.0 

0.0 


0.0 

0.0 


0.0 
0 . 0 


0 . c 

0 . G 


0.0 

3.0 


0.0 

c.o 


ALPHA 


c.c 

o.c 


l.oO 

0.0 


2.00 

0.3 


3.00 

0.0 


4.00 

0.0 


5.00 

0.0 


6.00 7 . 00 

0.0 0.0 


3.00 

0.0 


9.00 

0.0 


10.00 

3.0 
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Figure (7-4). Output data plot for Example 7-B 
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Table ( V— 5 ) . Output variaDlas of LAER0 1 



AL 


Angle of attack 


CLTOT 


Total coefficient of lift 


CDTOT 


Total coefficient of drag 


CLWP 

CLBW 

CLTP 


Wing panel coefficient of lift 
Additional lift on body due to wing 
Tail panel coefficient of lift 


CL3T 


Additional Lift on body due to tail 


CLB 

CDI 


Body alone lift coefficient 
Induced drag coefficient 


CNWP 

CNTP 


Wing panel normal force coefficient 
Tail panel normal force coefficient 


CLTD 


Lift coefficient due to tail deflection 


C DTD 


Drag coefficient due to tail deflection 


CN 


Total normal force coefficient 


CA 


Total axial force coefficient 


XCPW 


Wing center of pressure 


XCPT 

XCP 

CM 


Tail center of pressure 

Total missile center of pressure 

Total pitching moment about C.G. 


C DO WBT 


Zero lift drag coefficient of wing-body- 
tail combination 


CDMISC 


Miscellaneous zero lift drag coefficient 


C DOT 


Zero lift drag coefficient of tail 


C DO W 
C DOB 


Zero lift drag coefficient of wing 
Zero lift drag coefficient of body alone 


C DP ROT 


Drag coefficient of body protrusions 


CDINL 

CDAFT 


Drag coefficient of engine inlet 
Drag coefficient of boattail region 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



There are many topics which may be the subjects of fol- 
low on work contained within this thesis. Although the four 
programs have been installed on the IBM 370 computer system, 
these four alone do not fully satisfy the original goal of 
this work: Provide a computer supplement to the Tactical 

Missile Conceptual Design handbook. Numerous additional 
focal algorithms are utilized in the design handbook which 
deserve the attention of a programmer. Of immediate inter- 
est are the areas concerning radar or infrared guidance 
systems, baseline configuration modelling and weight distri- 
bution, and initial control and lifting surface design. 

Each of these topics can be programmed to provide missile 
design students interactive learning tools when coupled with 
the design handbook. 

The most urgent follow on work to this thesis is the 
restoration of the program LAER01 to a reliable, useful pro- 
gram. The program was modified and set up on the IBM 370 
computer system during the period immediately following the 
system's installation at the Naval Postgraduate School. As 
could be expected, the computer suffered many and varied 
growing pains in its early life. As a result of this, or of 
the human manipulation expanding the capability of the pro- 
gram, the effectiveness of LAER01 was substantially reduced. 

Work involving the other three programs would involve 
simply expanding their capabilities. The trajectory models 
program, LPATH, presently considers only two guidance laws: 
line-of-sight guidance and proportional navigation guidance. 
Other guidance laws which can easily be included in the pro- 
gram include pursuit, beam rider, and combinations of diffe- 
rent laws. It might also prove useful to be able to 
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simulate the entire missile trajectory but still only output 
the terminal phase of the encounter. Another option would be 
to provide the target with a controlled trajectory instead 
cf the constant acceleration condition now imposed. 

The propulsion sizing program, LPEOP, should have the 
various nozzle options incorporated into the program so that 
it isn't necessary to manually juggle the program output. 
Other booster-sustainer grain configurations could be 
explored, such as the booster grain being cast within the 
cere of the sustainer, or even a motor with only a single 
grain. Another suggestion for the convenience of the pro- 
gram users is to institute a shopping list of availlable 
propellants and their characteristics into the program. 

A similar list of available explosives and case materi- 
als could be put into the warhead sizing program, LBOMB. 
These shopping lists would provide ready access to reference 
information and, at the same time, decrease the number of 
data values to be manually input into the computer. Since 
the current program is limited to cylindrical warheads, an 
area cf expansion would be the flexibility of warhead 
styles, such as curved, shaped charge, continuous rod, etc. 

At present, the programs are somewhat hindered by the 
mechanics involved in producing the printouts and the plots. 
Due to the results of tailoring a program to be interactive, 
often it is required to completely exit a program before 
output can be received. Subsequent design iterations 
require re-entry into the program, which produces a certain 
justifyable annoyance to the user. Additionally, the compu- 
ter center has instituted a new policy of cancelling any 
jobs with duplicate job names, which can be severely irri- 
tating and cumbersome to the persons running the plot rout- 
ines contained within this thesis. However, the computer 
center has developed procedures which have the potential to 
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alleviate both of these problems. According to Volume 13, 
Number 4 of the Computer Center Newsletter . CMS commands can 
now be invoked from within a FORTRAN program. The print and 
plot operations presently contained within executive rout- 
ines can now be placed directly within the source programs. 
This will remove all current restrictions placed on- the num- 
bers of printouts received per session and will label each 
plot with the user's job name and not the project's job 
nase . 
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APPENDIX A. TRAJECTORY MODELS PROGRAM FLOWCHART 



MAIN PROGRAM ] 



Declarations 

DETORA (X) = ( X/3 60) *2*3. 14 16 
RATODE (X) =(X/(2*3. 1416) ) *360 

NRUN=0 

ISAVE=1 

I” I BAC K=0 1 



Input data: TIT LE. DEL, NAV , LO S 

IRT,TSP,TAC,TAL ,MSP,IRM 
■ 

r BACK=1 U 

:::c: 



f Call SCREEN 1 

::::cr 

r ipro=i i 

I LOS = 1 



Playback input data to terminal 
for verification 




TT3 = L0 S+TAL 
MTH = LOS+M AL 



LOS=DETORA| 


[LOS) 


TAL=DETORA < 


TAL 


TTH=DETO RA | 


TTH> 


MAL=OETORA 


MAL) 


MTH=DETORA | 


MTH) 



'f 
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R 1 = I RT 
P 1 =L OS 
A1=T AL 
ri=TTH 
R2=I RM 

:::r 



A2=ar sin (sin (A1 ) *TSP* IRM/ (!3SP*IRT) ) 
M a L ™ A2 
Mra=MAL+LOS 
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[ 



.1 



T3=TAC/TSP 
P3=TSP*sin (A 1) /R 1 



t 



] 



A3=T3 -P3 
R3=TSP*cos (A 1) 
R4=M3P*cos (A2) 






D6 = R4*rSP*sin ( A 1 ) 
D7=R2*TSP*A3*cos (A1) 
D8=-R3*MSP*sin (A2) 
D9=R1*MSP*COS (A2) 



A4= (D5+D7 + D8) /D9 
T4 = A4 + P3 



T 1=T1«-DEL*T3 
P1=P1 +DEL*P3 
A 1 = A 1 +DEL*A3 
R 1 = R 1 + 0 EL*R3 



i 



R2=R2+DEL*R4 
A 2= A 2 + 0 EL* A 4 
T2=r2+DEL*T4 



A (I) =H SP*T4 



[; 



i: 



T2 = A2 + P1 
R I) =R 1-R2 



: 



] 



XM (I) =XH (1-1) + D EL*M SP*cos (T2' 
YM (I) =YM (1-1 +DEL*MSP*sin (T2 
XT (I) =XT (1-1 + DEL*TSP*COS (T1 
YT (I) =YT (1-1) +DEL*TSP*sin (T1 

::::::: 

D (I) =D (1-1) +DEL 

1 
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f T5=DEL* (1-1) 1 
f Call SORT 1 
[ Call OUTPUT J 




A2=arsin (sin (A1) *TSP/msp) 
HAL- A2 
MT3 = MA L+LOS 



1 2 = M TH 

T 

[. 1 : 1 ] 





0 (I) 


RT-IRM 

i =0 
i =0 






1 


t 




km (: 
y m j: 

XT (I) =R (I 
IT (I) =R (I 


Fjio 

*cos (PI 
*sin (PI 


! 


i 


r 
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uo 



i 

R (I) =R (1-1) +D EL*R2 
A (I) =A (I- 1 ) +DEL*A3 
P 1 = P 1 +DEL*D2 
T 1 = T 1 *■ DEL*D6 
T2=T2+DEL*D7 
A 1 = A 1 <-0 EL*P2 
A2=A2+DEL*P3 



D (I) = D (1-1) * DEL J 




YES 




T5=DEL* (1-1) 

f Call SORT 1 

Call OUTPUT 1 

f Call SCREEN "U 

Do you want to run 
same problem again 



{yes 



the\> 

in? 



L0S=R ATODE (LOS) 
TAL=R AIODE (T AL 
MAL=RATODE (MAL) 



L 



1. 



< Do you want to rua\ 
a new problem? 

Jmo 

I" NR0N = 99 1 

T 

V.hU 



YES 
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APPENDIX 3. TRAJECTORY MODELS PROGRAM LISTING 



Following the next page is the program listing for the 
Trajectory Models program. It consists of three main divi- 
sions; the executive routines, the FORTRAN IV computational 
program, and the FORTRAN IV plotting program. The executive 
routines are used to establish required file definitions, 
initiate operation of the computational program, supervise 
the transfer of data to the plotting program, and provide 
operational information to the program user at appropriate 
times . 

The computational program, LPATH FORTRAN, consists of 
four subprogram divisions. The MAIN program accepts the 
input data and performs the calculations for the line of 
sight and the proportional navigation problems. Subroutine 
SCREEN is used to prompt the user to clear the terminal 
screen for proper positioning of the output. Subroutine 
SORT determines the largest missile acceleration value and 
the value ranges of the X and Y position coordinates for 
plotting reference. Finally, subroutine OUTPUT formats the 
calculated data and provides it to the user, the printer 
file, and the plot program data file. 

The plot program, PATHPLOT FORTRAN, is structured for 
direct submission to the MVS batch reader from the 3278 ter- 
minal. No cards need to be punched or read. Attention must 
be given to the correct JCL accounting data in the first 
line of the program; those shown in the listing are for 
illustrative purposes only. This program receives the data 
from the computational program and produces a single geo- 
graphic chart of the encounter in the encounter plane. Mul- 
tiple problems, up to nine, will overlay on the single 
chart . 
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FILE: LPATH EXEC A NAVAL POSTGRADUATE SCHOOL 

FILEDEF 08 DISK LPATH OUTPUT AO 
FILEDEF 07 DISK PATHPLOT DATA AO 
&BEGTYPE 

YOU WILL HAVE THE OPTION TO OBTAIN BOTH A HARDCOPY PRINTOUT 
AND A VERSATEC PLOT OF UP TO NINE ENCOUNTER GEOMETRIES. 

THE PLOT IS A SINGLE FRAME WITH ALL NINE GEOMETRIES 
SUPERIMPOSED ON ONE ANOTHER. THE HARDCOPY PRINTOUT IS 
IDENTICAL IN FORMAT TO THE TERMINAL OUTPUT. 

SEND 

LOAD LPATH 

START 

SBEGTYPE 

TO OBTAIN A HARDCOPY PRINTOUT OF THE RESULTS, TYPE AND 
ENTER: 

lpathpr 

THE OUTPUT WILL BE PRINTED ON THE VM PRINTER IN ROOM 140 
AND WILL BE IDENTIFIED B 1 YOUR USER NUMBER. IT USUALLY 
REQUIRES SEVERAL MINUTES TO OBTAIN THE PRINTOUT. 

TO OBTAIN THE PLOT OF YOUR ENCOUNTERS, TYPE AND ENTER: 

lpathpl 

THE PLOT WILL BE PRINTED IN THE COMPUTER ROOM AND WILL BE 
PLACED ON TOP OF THE ALPHABETIZED OUTPUT FILE IN ROOM 140. 
IT WILL BE IDENTIFIED BY THE JCL JOB NAME "PATHPLOT" AND 
USUALLY REQUIRES MANY MINUTES fEVEN HOURS!) TO BE 

OBTAINED. NOTE IF MANY USERS ARE REQUESTING PLOTS 

SIMULTANEOUSLY, THE COMPUTER CENTER PERSONNEL WILL CANCEL 
"EXCESS" JOBS USING THE AME IDENTIFIER. 

SEND 



FILE: LPATHPL EXEC A NAVAL POSTGRADUATE SCHOOL 

COPY LPATH PLOT A PATHPLOT DATA A PLOT FORTRAN A 
EXEC SUBMIT PLOT FORTRAN A 
ERASE PLOT FORTRAN A 



FILE: LPATHPR EXEC A NAVAL POSTGRADUATE SCHOOL 

PRINT LPATH OUTPUT (LINECOUN 70 
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APPENDIX C. WARHEAD DESIGN PROGRAM FLOWCHARr 



£ MAIN J 
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Subroutine 7 EL 



[ 



FRARE A= (FRMLB (J) /SDEN) **( 2 / 3 ) 
R= (TTH*HAr ) / (FRAREA**. 5) 



] 



BLV (J) +1 . 02546* (143 1 . 6 875*R) - (564. 1 857*R*R) 
+ (136. 7064*R**3) - (8. 774 4 7*R**%) 

BLV (J) = BL7 (J * (. 26/CDEN) 



K=RHO*FRARE A*0 . 65/FRMLB ( J) 
[I) * exp ( K*R KI LL) 
>J) /SOOND 



V(J) = VHIT 
M=V 
T=(1 + . 



!*M**2) *TAIR 




£ Subroutine SIZER "j 




f 
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- t_ 

FRAGS= (CDEN*5#PI/FRMLB (1 ) ) * (2 * RW*RW*CTH- R W*CTH*Cr H) 
FRLEN= (FRMLB ( 1 ) / (CTH*CDEN) ) ** . 5 
HIT3 = FR AGS* A 



£ Return 
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APPENDIX D 



WARHEAD DESIGN PROGRAM LISTING 



This program has two major sections; the executiv 
ine and the FORTRAN 17 computational program. The ex 
routine establishes the required file definitions and 
ates operation of the computational program. 

The computational program, LBOMB FORTRAN, consist 
four subprogram divisions. The MAIN program accepts 
input data, calculates the atmospheric characteristic 
formats and displays the output to the user and sends 
the printer file. subroutine 7EL calculates the init 
velocity required to propeL a given mass a specified 
tance through the atmosphere with a particular residu 
ocity remaining. It also determines the ballistic li 
velocities for the situation. Subroutine SIZER sizes 
warhead for a given Pd value. It also produces the c 
to-mass ratio, the number of fragments, the fragment 
and the average number of hits received by the target 
Subroutine SCREEN prompts the user to clear the termi 
screen for proper positioning of the displayed data. 



e rout- 
ecu + ive 
initi- 



s of 
the 

s , and 
it to 
ial 
dis- 
al vel- 
mit 
the 
harge- 
size, 

nal 



/ 
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FILE: LBOMB EXEC A NAVAL POSTGRADUATE SCHOOL 

FILEDEF 08 DISK LBOMB OUTPUT AO (PERM 
SBEGTYPE 

YOU MILL HAVE THE OPTION TO OBTAIN A HARDCOPY PRINTOUT 
OF THE LAST SOLUTION THAT YOU SOLVE. YOU MAY RERUN THE 
PROGRAM AS OFTEN AS YOU WISH BUT THE LAST RUN IS THE RUN 
THAT IS RECORDED. 

SEND 

LOAD LBOMB 

START 

&BEGTYPE 

TO OBTAIN A HARDCOPY PRINTOUT OF THE RESULTS. TYPE 
"PRINT LBOMB OUTPUT" AND ENTER. THE OUTPUT WILL BE 
PRINTED ON THE VM PRINTER IN ROOM 140 AND WILL BE 
IDENTIFIED 3Y YOUR USER NUMBER AND LABEL NAME. IT 
USUALLY REQUIRES SEVERAL MINUTES TO OBTAIN THE 
PRINTOUT. 

SEND 
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APPENDIX E 



PROPULSION SIZING PROGRAM FLOWCHART 



£ MAIN PROGRAM 
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i 

T=(PCB*DB^/J2* YIELD) 
APB=ATB/J 



[ 



I 



AM= PI* (DB/2-T) **2 
GEA=AM-APB 



] 



BRB = 1 .25 
B= (APB/PI) **. 5 
WT= (DB/2) - B 



ABB=DM/ (DENSB*BRB) 
RIB=BRB*TBB 




LB=VGB/GEA +DB/6+0 . 5 
WCB=2*PI DB/2 *LB*T*DENSC 
VOLB= AM* (LB-3 . 5) - ( . 0 025 *DB**4 ) 
NLB= (DEB-DrB) / (2*tan(ALN) ) 
TMB= HP B+WCB 




TBS= (R-R B) /VBF 
ALTS = ALT B+L ALT 

ELS=arctan ((TALT-ALTS) /(R-RB) ) 
ALTS=VBF*sin (ELS) *TBS*0. 5+ALTS 



f Call P RATIO 

I 
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i. 



BRS=0 . 9 9*BR S 



:] 



w= 


(2*»TS) /D 


B = e 


d (3. 509* W] 




C=( 


B-lt /(b+i; 




APS= ( 


1 . 065*0 1' 


/2 


GE 


A = AM - APS 



f r FAB 2=1 "L_ 

j EA= AM p— 

{ 



VGS=WPS /DENSS 1 

L l S=V"S/GEA J 




W T S = LS 

BRS= (AB2/AM) *BRS 
AB2= AM 



LS=LS+D/6*0.5 
WCS=2*PI*D/2*LS*T*DENSC 
VOLS = AM* (LS -0.5) -{. 0025 *D**4) 
TWS=WPS+WCS 
AES = AEB 
TS = T 




TBS= (R-RB) /VBF 
ALTS= ALTB+LALT 

ELS = ar ctan ( (r ALT- ALTS) / (R-RB) ) 



[ Call P RATIO 

f FS=DRAGS*1.U (WL-TWB) *sin (ELS) 

t 
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[ PRATIO ] 




£ PR= ( 1- (6. 8754e-6*ALT) ) **5.256 



I" PR = . 2234 *exp (-4. 80 6e-5* (ALT- 36 089) ) 

::::r 



■^4" P A=P R*1 4.696 1 

::x:: 

£ Return J 



132 



[ P3RAT ] 





B= 1- (PA/PC) 

C = 2* (G*G/ (G-1 )) * (2/ (3 



** ((3-1) /3) 

S + T) ) **( (G+1)/(G-1) ) 
.5} + ( (PE-PA 



CF = ,0.96.,(UcosMALSn/2,^B^*..| 

AI= (Cv*DM) / ( PC* 3 2. 174) 



* AR/P C) 



DT=2* (AT/3. 1415 9) **.5 
' * * **] (g-1) /G) -1) ) ** .5 



HE=( (2/ (G-1 ) ) * (PC/PE) **( (G-1) /G 
AE1=T+((3-1 )/2) *HE*KE 
AE= ( (1/MS) *((2/(G+1)( *AE1)**((3+1) 



D E=2* (AE/3 . 14159) 



n /(i 
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APPENDIX F 



PROPULSION SIZING PROGRAM LISTING 



Following this page is the program listing for the Pro- 
pulsion Sizing Program. It has two segments; the executive 
routine and the FORTRAN 17 computational program. The exe- 
cutive routine establishes the required file definitions and 
initiates operation of the computational program. 

The computational program, LPROP FORTRAN, consists of 
six subprogram divisions. The MAIN program accepts the 
input data from the user and perforins the guiding calcula- 
tions for the booster and sustainer motors. Subroutine 
PRATIO determines the ambient pressures at the design alti- 
tudes. Subroutine PCRAT defines the optimum chamber pres- 
sure to ambient pressure ratio with respect to the case 
material properties. Subroutine ARAT solves for the area 
ratio of the nozzle and tries to size the nozzle to the mis- 
sile diameter by varying the chamber pressure, characteris- 
tic velocity, and thrust coefficient. Subroutine SCREEN is 
used to prompt the user to clear the terminal screen for 
proper positioning of the output. And subroutine OUTPUT 
formats the calculated solutions and provides them to the 
user and to the printer file, if so directed by the user. 
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FILE: LPROP EXEC A NAVAL POSTGRADUATE SCHOOL 

FILEDEF 08 DISK LPROP OUTPUT AO (RECFM VA BLOCK 133 PERM 
&BEGTYPE 

YOU HILL HAVE THE OPTION TO OBTAIN A HARDCOPY PRINTOUT OF 
AS MANY ALTERNATIVES AS YOU WISH. THE PROGRAM WILL ASK 
YOU IF YOU DESIRE TO SAVE A PARTICULAR RUN, SIMPLY ANSWER 
ACCORDINGLY. 

SEND 

LOAD LPROP 

START 

&BEGTYPE 

TO OBTAIN A HARDCOPY PRINTOUT OF THE RESULTS, TYPE "PRINT 
LPROP OUTPUT" AND ENTER. THE OUTPUT WILL 3E PRINTED ON 
THE VM PRINTER IN ROOM 143 AND WILL BE IDENTIFIED 3Y YOUR 
USER NUMBER AND LABEL NAME. IT USUALLY REQUIRES SEVERAL 
MINUTES TO OBTAIN THE PRINTOUT. 

SEND 
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APPENDIX G . AERODYNAMIC COEFFICIENTS PROGRAM LISTING 



This program is divided into three major subdivisions; 
the executive routines, the FORTRAN II/IV computational pro- 
gram, and the FORTRAN IV pLotting program. The executive 
routines are used to establish the required file defini- 
tions, initiate operation of the computational program, 
supervise the transfer of data to the plotting program, and 
provide operational information to the user at appropriate 
times . 

The computational program, LAER01 FORTRAN, consists of 
five subprogram divisions. The MAIN program accepts the 
input data, conducts calc lations not done by other subrout- 
ines, formats the calculated data and provides the output 
data to the user, the printer file, and the plot program 
data file. Subroutine GE05UB calculates the missile wetted 
area and the Reynolds number per foot based on the flight 
altitude. Subroutine CLAS3B calculates the aerodynamic sur- 
face lift-curve slopes. Subroutine CATSUB calculates center 
of pressure locations, cross-flow drag coefficients, and 
interference factors. Subroutine SCREEN prompts the user to 
clear the terminal screen for proper positioning of the out- 
put. 

The plot program, AEROFLOT FORTRAN, is structured for 
direct submission to the MVS batch reader from the 3278 ter- 
minal. No cards need to be punched or read. However, 
attention must be given to the correct JCL accounting data 
in the first line of the program; those shown in the list- 
ing are for illustrative purposes only. This program 
receives the data from the computational program and pro- 
duces a group of six charts for each Mach number entered. 

The plots represent the relationships of Cl to alpha. Cm, 

Cn, Cd, Ca, and Cdi. The program can produce up to 24 sets 
of plots for a single run. 
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FILE: LAEROI EXEC A NAVAL POSTGRADUATE SCHOOL 

FILEDEF 08 DISK LAEROI OUTPUT AD (RECFM VA BLOCK 133 PERM 
FILEDEF 07 DISK LAEROI PLOT AO (RECFM VA BLOCK 80 PERM 
SBEGTYPE 

YOU HILL HAVE THE OPPORTUNITY TO OBTAIN BOTH A HARDCOPY 
PRINTOUT AND A SET OF PLOTS FOR ONE SET OF INPUT PARAMETERS 
EACH TIME YOU ENTER THIS PROGRAM. THE PROGRAM MAY BE 
RE-RUN CONTINUOUSLY AND YOU WILL HAVE THE OPTION TO CHANGE 
INPUT PARAMETERS FOR EACH SUCCESSIVE RUN, BUT YOU CAN 
OBTAIN THE PRINTOUT AND PLOTS PERTAINING TO THE LAST RUN 
ONLY. IF ADDITIONAL OUTPUT IS REQUIRED, RE-ENTER THE 
PROGRAM. 

SEND 

LOAD LAEROI 

START 

SBEGTYPE 

TO OBTAIN THE HARDCOPY PRINTOUT OF THE DATA TABLES, TYPE 
AND ENTER: 

LAEROI PR 



TO OBTAIN THE VERSATEC PLOT OF THE TABULAR DATA, TYPE AND 
ENTER: 

LAEROI PL 

SEND 



FILE: LAEROI PR EXEC A NAVAL POSTGRADUATE SCHOOL 

PRINT LAEROI OUTPUT (LINECOUN 70 
SBEGTYPE 

THE OUTPUT WILL BE PRINTED ON THE VM PRINTER IN ROOM 140 
AND WILL BE IDENTIFIED WITH YOUR USER NUMBER AND LAST NAME. 
IT USUALLY REQUIRES SEVERAL MINUTES TO OBTAIN THE PRINTOUT. 
SEND 



FILE: LAER01PL EXEC A NAVAL POSTGRADUATE SCHOOL 

COPY LAEROI PLOT A LAERO PLOTDATA A PLOT LAEROI A 
EXEC SUBMIT PLOT LAEROI A 
ERASE PLOT LAEROI A 
SBEGTYPE 

THE PLOT WILL BE DRAWN IN THE COMPUTER ROOM AND PLACED OVER 
THE ALPHABETIZED OUTPUT FILE IN ROOM 140. IT WILL BE 
IDENTIFIED BY "AEROFILE" AND USUALLY REQUIRES MANY MINUTES 
(HOURS, DAYS) TO OBTAIN. 

SEND 
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